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ABSTRACT 


A quasl-geostrophlc numerical model, including friction, 
radiation, and the observed planetary topography, is applied 
to the general circulation of the Martian atmosphere in the 
southern hemisphere at latitudes south of about 35° . Near 
equilibrium weather systems developed after about 5 model days. 

To avoid violating the quasi-geostrophic approximation, only 
0.8 of the already smoothed relief was employed. V/esther 
systems and velocity fields are strikingly tied to topography. 

A 2mb middle latitude jet stream is found of remarkably terrestrial 
aspect. Highest surface velocities, both horizontal and vertical, 
are predicted in western Hellas Planltla and eastern Argyre 
Planltia, which are observed to be preferred sites of origin 
of major Martian dust storms. Mean horizontal velocities > 30 m s” 
and vertical velocities > 0.2 m s”^ are found just above the 
surface velocity boundary layer. When consideration is taken of 
scaling to full topography and the probable gustiness of Martian 
winds, it seems very likely that the general circulation is 
adequate, at certain times and places, to transport dust from 
the surface of Mars, as observed. Certain sources and sinks of 
vertical dust streaming are suggested; the entire south cir- 
cumpolar zone appears to be a dust sink in winter. 


1. Introduction 


Leovy and Mintz (1969) constructed the first numerical model 
of atmospheric circulation on Mars. They employed a primitive 
equation two-level model with radiative-convective heating and 
cooling. Allowance was made for the condensation of carbon dioxide 
on the surface and the consequent release of latent heat. At 
winter solstice, maximum meridional winds of about 8 m/s for bo*-h 
atmospheric levels were found; at equinox, the corresponding 
values were about 2 m/s. On the other hand, the maximum zonal 
winds at high p^tltude at solstice reached about 70 m/s at ^5® 

North latitude, while at equinox the maximum zonal winds were 
about 25 m/s. The corresponding maximum "surface" winds — 
corresponding to the region just above the surface boundary 
layer — were estimated to be about 20 m/s at solstice and about 
10 m/s at equinox. No topography was included. 

However, the Martian surface is know to exhibit major elevation 
differences (Sagan et al., 1967; Sagan and Pollack, 1968; Pettenglll 
et al., 1973; Downs et al., 1973; Conrath et al., 1973; Kllore 
et al., 1973; Hord et al., 197^). Over continental dimensions, 
elevation differences comparable to a scale height occur on Mars; 
and over scales of a few hundred kms elevation differences of 2 or 
3 scale heights exist. Under these conditions, winds produced by 
relief seem able in restricted locales on Mars to Increase the mean 
winds by factors of two to three over the mean wind velocities 
expected in the absence of relief (Gierasch and Sagan, 1971). 
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Wind velocities at half surface pressure level 
in excess of 8o ra/s were therefore predicted at certain seasons 
and locales, and a range of observational consequences of such 
high winds predicted (Sagan et al., 1971). Time-variable streaks 
and splotches on the Martian surface have been interpreted in 
terms of windblown dust produced by such high winds (Sagan, et 
al., 1972; 1973; 197^; Veverka et al., 1974). 

The work reported in the present paper is a simple numerical 
experiment on the general circulation of the Martian atmosphere 
with topography included. The analysis is restricted to the 
Southern hemisphere only. The assumed topography is shown in 
Figure 1, and is based upon Mariner 9 UVS, IRIS and S-Band 
occultation experiments. Near the South Pole the assumed eleva- 
tions reflect some Mariner 9 television stereoscopic results as 
well (Blaslus, 1974). Since some inconsistencies in the results 
of the various Mariner 9 and groundbased radar topographies still 
exist, the Southern hemisphere Martian topography adopted in 
these calculations may not correspond precisely to the true cir- 
cumstances. Nevertheless, the general relief seems to be 
sufficiently well-established that the Influence of topography 
on atmospheric circulation can be tested. 
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2. Structure and Equations 

The physics of this two-level model is based upon the 
assumption of a hydrostatic, adiabatic, quasl-geostrophlc, ideal 
gas atmosphere entirely composed of carbon dioxide. Radiative 
transfer between the atmosphere, space and ground is Included 
in the model, as is surface friction. No provision is made 
for the latent heat of phase change of atmospheric constitutents 
or for albedo differences over the Martian surface. 

The model's vertical structure, shown in Table 1, is one 
in which the geopotential heights of two pressure surfaces, at 
2 and 4 mb, are predicted, while at the 3mb level the vertical 
velocity is diagnosed. At the 6mb surface an orographic 

vertical velocity is computed using the surface topographic 

gradient and the 4mb surface horizontal velocity. Finally, at 
the upper surface the vertical velocity is assumed to be zero. 
The set of four predictive and diagnostic equations used in 
this model can be derived using the hydrostatic, adiabatic, 
quasi-geostrophlc and ideal gas approximations (see, e.g., 
Haltiner, 1971; Thompson, 1961). 

First, at the upper, 2rab, surface the partial differential 
equation for predicting changes in the geopotential heights 


of that surface is 



( 1 ) 
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where J is the Jacobean operator; f Is the Coriolis parameter 
■ 2nsin6; 6 is the latitude; d is the distance between grid 
points; m is the map scale factor which corrects the distances 
for the projection used (in this case polar stereographic) 
to the real world values; is the geostrophic relative 

vortlcity on the 2mb surface ■ (g/f) V Z2mb* difference 

in pressure between the two geopotential surfaces (2rab); and 
W3mb vertical velocity of the 3mb surface in units of 

pressure/time. 

On the lower 4mb surface the equation for geopotential 
change is similar but has additional terms including friction 
and the surface orographic vertical velocity, 
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Here g is the gravitational acceleration » 372 cm s and P 

y 

and are the y and x frictional components at the surface. 
They are treated in a manner similar to that employed by Leovy 
and Mintz (1969), l.e., 
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where (p^-p*.) Is the difference In pressure between the upper 

9 V 

and lower levels. The surface stress is written 

Is - <"> 

where the friction velocity U* Is 


U* * 



(5) 


Y Is the horiijntal wind at the surface (6mb), which Is derived 

'wg 

from a linear extrapolation of the winds at the 2 and 4mb levels, 
2 

and C Is the momentum drag coefficient. Assuming a stable 
m 

stratification, it has the value 0.9 x 10”^ (cf. Leovy and 
Mints, 1969). p. is the atmospheric density at this surface. 

This lower boundary friction is the principal means of 
dissipation of thermal energy. 

The vertical velocity of the 3mb surface is diagnosed 
and the radiative effects are introduced into the model by 
means of the omega equation: 


— , 2 

“ Sf^IT' ^ ^2mb* ^^4mb ” '^2mb^ 

“W '^•^^^Umb* ♦iimb " •^^^2mb* ♦2mb 
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o is the static stability ■ (gp6}~^ (36/dp) ; 6 Is the potential 
temperature ■ TCp^/p)*^, where < ■ Cp/R; P^/P ■ (6mb/3mb) • 2; R is the 
gas constant and Cp Is the specific heat of CO^ at constant pressure . 

A value of -^1. 1^32*1 x 10“^ m^ erg"^ was employed for o, 0.188 
X 10^ ergs gm"^ for R and 0.8?^ x 10^ ergs gm“^ (®K)“^ for Cp. 

Q, the heating or cooling rate is written as a simplified 
version of the radiative term found in Blumsack, Qlerasch and 
Wessel's (1973) model of the Martian atmospheric structure: 


® * ^solar * '^planetary 


(7) 


The solar heating is found from the relation 


- 0.216 X 10^ 65.30 ®K cm"^sec“^. 


where arccos C is the solar zenith angle. 
The planetary term Qp is 


q_ - (§-) C JT ][t— Vv— ] » 

p p (1-s^T + D 

C PS /d ^ 


( 8 ) 


Here, a • Hi the two radiation parameters are C ■ 55 cm 
—2 

and D • 383 dyne cm ; and B is the Planck function at the 
center of the 15u band, i.e.. 
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® t"® cn'^seo"^] (9) 

Bn, BCO'*') and B(l~) denote the Planck function of the ground 
and the top and base of the atmosphere, respectively. Thus, the 
first term of eq. (8) describes the heating or cooling by the 
ground and the second term describes heat loss into space. 

At the lower boundary the orographic vertical velocity 
is computed from 


"6mb ■ - 

7h is the gradient of the topography and is the horizontal 

velocity on the 4mb surface. 

Finally, in order to diagnose the mid-atmospheric temperature 
throughout the simulation the thickness equation was used, i.e.. 


where T is the average temperature between the 2mb and 4mb 
surfaces and Az is the altitude difference between the two 


surfaces. 
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3* Numerical Methods 


The simulations were made on a 33 x 27 grid covering the 
majotlty of the non-equatorlal area of the southern hemisphere 
of Mars. The grid spacing was 202.85 km and the projection used 
was a polar stereographic, true at 60®S. The omega equation and 
the two geopotential predictive equations were solved by a 
’^.lebmann overrelaxation method In which the computer continuously 
guesses at a solution, tests If the error Is below a certain 
tolerance, and, if it Is not, subtracts or adds by an amount 
slightly greater than the actual error (cf. Haltiner, 1971; 
Thompson, 1961). 

For stability a time step of 30 minutes and an Arakawa 
Jacobean was used. Dealing with the boundaries proved more 
difficult but excellent results were realized by using cyclic 
boundaries In which surface features could pass off one edge 
to reappear on another. This was effected by first splitting 
both the outside rows of the active simulation area and those 
TDws adjacent and exterior to them (see Figure 2). Then the 
following procedure was adopted for each of the four sides: 
each outside half row (not in the active simulation area) 
takes on the values of the diagonally opposite half row located 
in the active simulation row adjacent to it. Thus in Figure 2 
point A takes on the parameter value of B, point C takes on the 
value at D, point E takes on the value of F, point M 
takes on the value at N, and so on through that side and the 
other three sides. To avoid boundary errors, all results 
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are restricted to latitudes poleward of 30°S, and the outermost 
latitude circle shown on our polar sterographic projection is 
40®S. 

The solar zenith angles at the various grid points were 
allowed to vary diurnally so as to more closely simulate the 
tl^ue heating. Furthermore, these angles could also be varied 
to simulate the seasons. All simulations begin with an isothermal 
atmosphere (winter, l85*^K; spring, 190°K; sunoner, 200®K), and 
with a latltudlnally prescribed ground temperature. 

The simulations were generally performed for 7*5 or 8.5 
days, the limiting factor being computing costs rather than 
model stability. On the Cornell IBM 360/65 an 8.5 day simulation 
took approximately 40 minutes of central processing time. 
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ll. Validity of the Results 

Since the frictional and radiative time constants of the 
Martian atmosphere are short, >2-3 days. It Is expected that 
a simulation of 7*5 or 8.5 days should give definitive results. 

It was In fact observed that In all simulations a near-equlllbrlum 
situation did establish Itself within 4.5 days Into the simulation 
(see Figure 3). 

Another question deals with the validity of the non-geostrophlc 
terms for heating, friction and, most significant, orographic 
vertical veloclf.y. The first two terms are not of large 
magnitude and should pose no problems. However, the effect of 
the large vertical orographic velocities on the validity of the 
model must be examined. We define a nondlmenslonal vertical 

velocity w* -2 w/W, where W is the vertical velocity scale - 

2 

fjjU /g. Here f^ Is the mean Coriolis parameter over the latitudes 
In question and U Is the mean horizontal velocity scale. When 

w* « Ro = U/f^L (12) 

where L Is a characteristic distance scale, the quasl-geostrophlc 
approximation remains valid. More specifically. It can be shown 
(Phillips, 1963) for orographic vertical velocities, that an 
additional condition for the validity of the approximation Is 
the following relation connecting the Rossby number, Ro, the 
atmospheric scale height, H, and the height, h^^, of the largest 
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topographlc feature: 

< CRo) H (13) 

max 

We adopt H “ 8 km. For Mars Ro = 0.25 - 0.50, although for 
large slopes (> 1 percent) the Rossby number approaches unity 
(cf. Glerasch and Sagan, 1973). These conditions are satisfied 
if we can assure smaller slopes with heights ranging less than 
one-half a scale height. Since the actual Martian topography 
is so rough that these conditions are not fulfilled, numerical 
methods of smoothing and scaling were employed. For most of the 
runs a linear relief scaling of 0.8 was used with a nine point 
smoother. This proved quite effective. As a test, however, we 
made a few runs without smoothing or scaling and found the model 
to be dramatically less stable; clearly the quasigeostrophic 
approximation was failing under conditions of realistic Martian 
relief. 

Because the quasigeostrophic approximation fails at 
equatorial latitudes, the calculation of heating rate does 
not include the solar radiation Intercepted by the Martian 
surface at low southern latitudes. To perform a correction, we 
assume that there is negligible Interhemlspherlc heat transfer 
and that there are no systematic dependences of albedo on 
latitude. (The latter assumption is in fact strongly contradicted 
by the data, and the region between 35® S and the equator has 
a marked concentration of low albedo regions). A straightforward 
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geometrlc calculation shows that the ratio of the insolation 
between 0 and 35®S to that between 35 and 90°S at winter noon 
is about 2,4. The meridional heat flux is, to first order, 
proportional to the product of the meridional wind and the 
equator-to-pole temperature gradient, and therefore to the 
square of the meridional wind. Thus, if the wind is proportional 
to the square root of the heat flux, our horizontal and vertical 
winds should both be corrected by a factor of 2.4^ a 1.5* 

Because of the topography truncation and the albedo assumption, 
the mean winds derived are probably lower limits on their 


actual values. 


5. Results 


The principal topographic features of Figure 1 are 
the prominent basin in Hellas Planltla (near 290^W, 59^3), 
a shallower elongated depression centered in Argyre Planitia 
(near 40®W, 45®S), and a depression offset from the pole (near 
110®W, 75®S). Hellas and Argyre are frequently reported to be 
connected with Martian dust storms, and one objective of this 
study was to examine topography-connected vfeather systems and 
velocity fields in these locales. Because of the possibility 
that there is a permanent dust pall in these basins (Sagan et al., 
1971) they may be substantially deeper than the depths adopted 
from IRIS and UVS Mariner 9 data in Figure 1, in which case 
the weather systems and velocity fields generated by these basins 
should be yet more pronounced. Figure 1 also displays extended 
plateaux on the Noachis/Mare Serpentls boundary (near 3^5®W, 

35® S) and in a broad region centered on Thaumasla Fossae, south 
of Solis Lacus (near 90®W, 40®S). While future refinements 
of the Martian topography may alter some of the details of our 
results, the general character of the results should remain 
valid. 

In Figure 4 are exhibited numerical results on the develop- 
ment of Martian weather systems in winter at the 4mb level. 
Geopotential heights are given for seven steps in the development 
of the system, from 0.5 days to 7.5 days, at which time the 


convergence of the calculations is evident. Weather systems 
whose positicxis ai^ar to be connected to t(^graphy in Figure % Include 
the high near Noachis, and the extension of the 3.6 km geopotential 
surface into Hellas. The weather system is also characterized 
by a stable polar low. The corresponding 7*5 day result for 
summer is shown in Figure 5a, where the high and both lows 
appear to have no direct connection with the topography. The 
corresponding 8.5 day result for spring equinox is displayed 
in Figure 5b. Here there is a pronounced low in Western Hellas; 
a high in Eastern Argyre; a high on the Noachls/Mare Serpent is 
plateau; and a low in the polar depression. 

In general, at the ^mb surface a regime of closed-off 
regions of low and high geopotential develop. Many of these 
weather systems, especially in the winter and spring simulations, 
are approximately stationary. That their positions were deter- 
mined by topography was conclusively established by varying the 
relief. A marked strengthening of stationary system was observed 
as the topography scaling factor was increased. These results 
should be compared with the topography-independent weather maps 
of Leovy and Mlntz (1969, Figures 19a and b). 

The 2mb weather maps of our calculations are shown for 
winter and summer solstices and for spring equinox in Figure 6. 

The corresponding horizontal velocities for the winter solstice 
map are exhibited in Figure 7. An upper level Jet stream with 
maximum velocities of H5 m s”^ is evident. Its location at 
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mlddle latitudes and the order of magnitude of Its velocities 
are both strikingly similar to the terrestrial case. Also 
analogous to the Earth's Jet stream Is the greater Intensity 
of the winds during winter, with wind maxima found In closed- 
off jet cores. 

Horizontal velocities at the level are shown In Figure 
8 at Spring Equinox for the cases of no topogr c ly and 0.8 the 
full topography of Figure 1. It Is apparent that Inclusion of 
topography greatly complicates the velocity field as well as 
Increasing the number of regions experiencing the highest 
velocities. This trend continues for surface velocities. 

In Figure 9 we see the Winter Solstice velocity field for 0.8 
full topography. The highest velocities shown are about 30 m s“^, 
and correlate very well with the locales of highest velocity 
for the 4mb surface (Figure 8b), In the spring. This correlation 
is another demonstration of the high degree of topographic control 
of Martian winds. Since we know that there are locales many 
km higher and lower than the smoothed topography of Figure 1, 
occasional velocities two or three times larger than those exhibited 
in Figure 9 may be Inplied. Of the five locales In Figure 9 of hi^st velocity, 

all but one are in the western Interior slopes of Hellas Planltla, 
encompassing the Hellespont! Montes; in the center of Noachis; 
and in the southeastern quadrant of Argyre Planltla. (The fifth 
Is on the border of Eridanla and Mare Tyrrhenum. ) But these 
are Just the regions in which large dust storms have been 
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observed over the course of a century of groundbased observations 
of the planet. While the calculations of Figure 9 are for Winter 
Solstice the strong degree of topographic control of the highest 
winds suggests that these will be the locales of highest winds 
for all seasons. Recent calculations of the horizontal winds 
above the velocity bou? xayer necessary to lift dust from 
the Martian surface range from about 25 to about 80 m s"^ 

(Sagan and Veverka, 1976; see also Sagan et al., 1973). The 
results of the calculations of the present paper appear to be 
In the required range, even without Invoking the tail of the 
velocity time distribution function (gusts). To the extent 
that the higher velocities mentioned are required for dust 
mobility, topographic generation of the highest winds Is even 
more strongly Indicated. 

Huss (197^) has summarized a range of available Information 
on terrestrial wind velocity distribution functions. Because of 
the steeper temperature gradients on Mars, terrestrial gustiness 
should provide a lower limit on Martian gustiness. The results 
imply that for a 15 m/s mean horizontal velocity on Mars there 

_-3 

will be a 10 probability of a gust, for a minute's duration 
and a mile of surface, of velocity 40 m/s. Extrapolation suggests 
that where the mean wind is 30 m/s there is a 10 probability of 
such a gust of velocity 70 m/s. Thus, a reasonable allowance for 
our topography truncation, forced on the calculations by the 
quaslgeostrophic approximation, and likely levels of Martian 
gustiness, imply that Martian general circulation topographic 
winds are entirely adequate to produce saltation and suspension 


of aerosols on Mars. 

The vertical velocities at the surface for winter, spring 
and sununer are shown in Figure 10. The vertical velocities 
for the same three seasons at the 3mb level are shown in 
Figure 11. The highest upward velocities in spring and summer 
are found in western Hellas/Hellespontl Montes, which is also 
a region of high upward velocity in winter. The highest winter 
velocities are found in northeastern Argyre. These highest 
upward velocities are about 0.1 to 0.2 m s“^. From a universal 
scaling between drag coefficient and Reynolds number, Sagan 
and Bagnold (1975) have calculated the modified Stokes fall 
velocities under a range of Martian conditions. The vertical 
velocities derived in the present paper then correspond to 
suspension of particles 10 to 20 pm in diameter. Such particles, 
falling out from an 
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altltude of 10 km, will reach the ground in about a day. Thus, 
their continued presence requires continued lifting by the winds. 

On the other hand, particles about 1pm In diameter have a 

7 

characteristic fallout time from 10 km altitude of about 10 s, 
which Is the typical duration of a major Martian dust storm, such 
as the one experienced by Mariner 9. Likewise 1pm Is the 
typical grain diameter deduced for the 1971 dust storm from 
a range of different techniques (see, e.g., Hartmann and Price, 197^). 
Thus it appears that the vertical and horizontal wind ccoipoients in Western Hellas 
and Eastern Argyre produce copious amounts of suspendible dust each year. The 
large particles fall out rapidly, not far from their sites or 
origin, but the smaller particles are carried over much of the 
planet before they sediment out. In general the vertical 
velocities on Mar? appear often to be high enough to cause 
suspension of pm diameter particles; the reason that the planet 
is not perpetually dust-covered is that the horizontal winds are 
only infrequently high enough to produce saltation and suspension. 
These calculations suggest that major and minor dust storms 
both should be a common occurrence on Mars, a conclusion in 
apparent accord with recent observational experience. 

Locales of preferred downward velocity are also of some 
Interest. For example, In all seasons eastern Argyre appears 
to be a source of suspendible dust, and western Argyre a sink. 

There may of course be interior circulation of dust within the 
Argyre basin (Sagan et al., 1971). The polar depression appears 
to be a source of suspendible dust in spring, and a sink in 
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suRimer and winter. The entire south circumpolar zone appears 
to be a dust sink In winter (Figure 10a) » a fact which may be 
related to the probably deposltlonal polar laminae (Murray et 
al., 1973). 
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TABLE 1 

Structure of the Model 
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Plgure Captions 

Figure 1: Polar stereographic projection of adopted Martian 

southern hemisphere topography. The outermost clrciunpolar 
circle corresponds to latitude 40^S; the circles here and In 
subsequent such diagrams are spaced In 10** latitude Intervals. 
The topography exhibited here, compiled from a variety of 
Mariner 9 data, Is smoothed, and considerably greater 
elevation differences actually exist. 

Figure 2: Parameterization scheme for cyclic boundaries In 

the active simulation area. See text for details. 

Figure 3: Stability of numerical solutions after several model 

days for three different Martian seasons. 

Figure Winter solstice weather systems at the 4mb level. 

The contours are geopotential heights in kms. Development 
of systems Is shown for elapsed times between 0.5 days (4a) 
and 7-5 days (4g). Highs and lows are Indicated. 

Figure 5: Weather systems at the 4mb level for summer solstice 


after 7.5 days (5a) and for spring equinox after 8.5 days 
C5b). 
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Plgure 6: Weather systems at the 2mb level for winter solstice 

after 7.5 days (6a), for spring equinox after 8.5 days (6b), 
and for summer solstice after 7*5 days (6c). 

Figure 7: Horizontal velocity field in m/s at 2mb level for 

winter solstice after 7*5 days. 

Figure 8: Horizontal velocity field in m/s at Mmb level for 

spring equinox after 7.5 days. Calculations are shown 
for no topography (8a) and for 80 J{ of the full (smoothed) 
topography adopted in Figure 1 (8b). 

Figure 9: Horizontal velocity field in m/s Just above velocity 

boundary layer for the Martian surface. Calculations 
are for winter solstice after 7.5 days and assume 80? of 
the full topography. 

Figure 10: Vertical velocity field in cm/s Just above velocity 

boundary layer for the Martian surface. Calculations are 
after 5.5 days for winter solstice (10a), spring equinox 
(10b), and summer solstice (10c). Solid contours indicate 
upward flow, dashed contours downward flow. 

Figure 11; Vertical velocity field in cm/s at the 3mb surface 

after 5.5 days for winter solstice (11a), spring equinox (lib), 
and summer solstice (11c). Solid contours indicate upward 
flow, dashed contours downward flow. 
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